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SURVEY 

At  the  St.  Louis  meeting  of  the  American  Physical  Society,  on  March 
25-28,  we  are  presenting  a  review  of  recent  work  on  the  energy  bands  of  solids, 
using  the  augmented  plane  wave  method.  I  am  giving  a  half-hour  invited  paper, 
followed  by  contributed  papers  by  Mattheiss,  Wood,  Switendick,  and  Pratt,  all 
on  topics  which  have  been  recently  reported  in  these  Quarterly  Progress  Re¬ 
ports.  There  are  several  points  which  I  shall  emphasize  in  my  talk,  as  well 
as  sketching  the  main  ideas  of  the  APW  method.  One  of  these  points  is  that 
the  APW  method  is  not  tied  to  the  use  of  a  potential  which  is  constant  between 
the  atomic  spheres,  as  well  as  being  spherically  symmetrical  within  the  spheres 
There  is  no  complication  involved  in  having  an  arbitrary  potential  in  the  space 
between  the  spheres.  All  that  has  to  be  done  is  to  expand  in  three-dimensional 
Fourier  series  a  function  equal  to  the  assumed  potential  between  the  spheres, 
and  zero  within  each  sphere.  The  Fourier  components  of  this  expansion 
immediately  give  contributions  to  the  diagonal  and  non-diagonal  matrix  elements 
of  the  Hamiltonian  between  the  augmented  plane  waves,  before  solving  the 
secular  problem,  which,  carries  through  as  usual.  This  was  pointed  out  in  the 
first  paper  on  the  APW  method,  and  has  been  emphasized  particularly  in  still 
unpublished  work  by  P.  M.  Marcus  and  H.  Schlosser  (I  am  indebted  to  them 
for  pointing  this  out  to  me).  Though  we  have  not  so  far  used  such  a  Fourier 
expansion,  we  are  actively  exploring  its  use,  and  shall  expect  to  use  it  as  soon 
as  we  have  decided  what  potential  to  assume. 

A  second  point  which  I  shall  emphasize  is  that  though  the  energy  bands 
in  some  metals  are  rather  free-electron  like,  nevertheless  this  is  not  a  general 
situation,  and  it  makes  the  use  of  the  various  pseudo-potential  methods  which 
have  been  proposed  very  dubious  in  general.  One  can  approach  the  problem  by 
considering  the  logarithmic  derivatives  of  the  radial  solutions  of  Schrodinger’s 
equation  at  the  radius  of  the  atomic  spheres.  If  these  logarithmic  derivatives 
happened  to  equal  the  logarithmic  derivatives  of  the  spherical  Bessel  functions 
corresponding  to  free  electrons  whose  energy  equalled  the  actual  energy  of  the 
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problem,  the  solution  of  Schrodinger's  equation  inside  the  sphere  would  join 
smoothly  onto  a  single  plane  wave  outside  the  sphere.  This  would  mean  that  a 
plane  wave  solution  in  the  regions  between  the  spheres  would  not  be  scattered 
by  the  atoms.  In  such  a  case  one  can  show  that  the  solution,  for  this  energy, 
is  free-electron  like  as  far  as  its  energy  is  concerned,  as  well  as  having  a 
wave  function  which  is  a  plane  wave  except  within  the  spheres.  This  is  the  same 
situation  which,  in  the  problem  of  scattering  of  a  plane  wave  by  an  isolated  atom, 
leads  to  the  Ramsauer  effect,  according  to  which  certain  atoms,  particularly  the 
inert  gas  atoms,  show  practically  a  zero  scattering  cross  section  for  slow 
electrons.  The  point  now  is  that  the  logarithmic  derivatives  for  the  crystal 
potentials  in  such  cases  as  sodium,  magnesium,  and  aluminum  have  almost 
exactly  this  property,  leading  to  the  free -electron-like  behavior  of  the  energy 
bands  in  these  cases.  In  most  cases,  however,  we  do  not  have  this  situation, 
and  the  use  of  free  electron  solutions  as  a  starting  point  is  not  at  all  justified. 

In  every  case  the  APW  method  is  applicable,  but  the  only  way  to  test  the  appli¬ 
cability  of  pseudo-potential  methods  is  to  compare  them  with  more  accurate 
solutions,  such  as  the  APW,  and  if  the  accurate  solution  is  available,  there  is 
no  point  in  using  the  free-electron  approximation. 

A  further  point  regarding  the  APW  method  is  its  applicability  in  many 
cases  where  other  methods  are  impossible  or  difficult.  One  of  these,  of  course, 
is  the  study  of  the  transition  series  of  elements.  Another  is  the  study  of  crystals 
involving  very  heavy  atoms,  in  which  no  particular  difficulty  arises  when  using 
the  APW  method,  as  Pratt  and  his  students  have  shown,  but  for  which  the  OPW 
method  is  rather  difficult.  For  the  study  of  spin-orbit  interaction,  the  APW 
method  is  particularly  convenient,  since  the  spin-orbit  interaction  Hamiltonian 
is  zero  except  where  there  is  an  external  electric  field,  which  means  that  if  we 
have  the  conventional  potential  used  in  the  APW  method,  constant  between  atoms, 
the  Hamiltonian  will  vanish  except  within  the  spheres.  Within  these  spheres, 
since  we  have  spherical  symmetry,  we  can  carry  through  the  ordinary  transfor¬ 
mation  reducing  the  spin-orbit  interaction  to  an  operator  involving  S*  L,  the 
scalar  product  of  spin  and  orbital  angular  momentum.  Furthermore,  the  wave 
function  within  the  spheres  is  expanded  in  series  of  terms  each  with  a  fixed 
angular  momentum,  so  that  it  is  easy  to  deal  with  the  matrix  elements.  These 
advantages  have  been  pointed  out  in  an  unpublished  memorandum  by  Larry 
Johnson,  working  with  Pratt. 

I  have  included  this  sketch  of  my  remarks  at  the  meeting,  since  they 
will  not  be  published  elsewhere.  As  for  the  new  results  to  be  reported,  per¬ 
haps  the  most  striking  are  Mattheiss's  calculations  of  energy  bands  for  almost 
the  whole  3d  transition  series  of  elements,  which  he  reports  in  the  present 
Quarterly  Progress  Report.  These  are  still  preliminary  results,  in  that  he  has 
so  far  determined  energy  bands  only  in  a  few  directions,  but  the  programs  are 
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working  so  well  that  we  can  now  turn  out  these  energy  bands  very  rapidly,  and 
it  will  not  be  long  before  much  more  complete  information  is  available.  It  is 
particularly  interesting  to  check  the  validity  of  the  rigid  band  model  of  the 
transition  elements,  to  see  the  gradual  narrowing  of  the  3d  bands  as  we  go  up 
in  the  series,  and  to  see  the  rapid  descent  of  the  3d  bands  below  the  4s  and  4p 
as  we  go  from  copper  to  zinc  and  gallium.  This  is  the  first  time  that  it  has  been 
possible  to  assemble  on  a  common  basis  a  set  of  energy  band  calculations  for 
these  elements,  including  representatives  of  body-  and  face -centered  cubic  and 
hexagonal  structures.  Wood's  work  on  gallium  fits  onto  this  series,  and  is 
entirely  consistent  with  it.  The  case  of  gallium,  however,  is  one  of  the  most 
complicated  which  we  have  yet  encountered  as  far  as  the  structure  of  the  Fermi 
surface  is  concerned,  and  it  is  by  now  clear  that  many  further  points  in  the 
Brillouin  zone  will  have  to  be  calculated  before  the  shape  of  the  Fermi  surface 
is  entirely  clarified.  Wood  is  proceeding  with  the  calculation  of  these  further 
points. 

In  the  two  preceding  Quarterly  Progress  Reports  I  have  presented  re¬ 
sults  on  symmetry  properties  of  a  number  of  hexagonal  and  cubic  space  groups. 

I  conclude  this  series  with  some  tetragonal  and  rhombohedral  space  groups,  in 
the  present  Quarterly  Progress  Report.  In  the  next  Ouarterly  Progress  Re¬ 
port,  of  July  15,  19^>3,  I  propose  to  give  some  discussion  of  double  space  groups, 
and  their  relation  to  spin-orbit  interaction. 

The  work  on  molecular  calculations  and  molecular  integrals  has  been 
proceeding  rapidly  during  the  preceding  quarter,  as  is  shown  by  no  less  than 
seven  contributions  to  this  Quarterly  Progress  Report,  by  Moskowitz,  Sutcliffe, 
Harrison,  Ellis,  Wright,  Dahl,  and  Nielson  and  Woznick.  Some  of  these  deal 
with  specific  molecules,  but  most,  of  them  are  concerned  with  the  calculation 
of  three-  and  four-center  integrals,  with  mathematical  techniques  associated 
with  their  calculation,  and  with  packages  of  programs  for  the  calculation  of 
molecular  wave  functions.  Excellent  progress  is  being  made  with  this  pro¬ 
gramming,  but  there  is  still  quite  a  way  to  go  before  we  are  ready  to  make 
wholesale  calculations  of  polyatomic  molecules. 

We  have  one  addition  to  the  group  since  the  last  Quarterly  Progress 
Report,  Dr.  J.  P.  Dahl,  from  the  University  of  Copenhagen,  who  will  be  with 
the  group  for  approximately  a  year.  Fred  Quelle,  who  has  been  associated 
with  the  group  for  several  years,  as  a  graduate  student  at  Harvard  as  well  as 
being  associated  with  the  Lincoln  Laboratory,  has  now  completed  his  require¬ 
ments  for  the  doctor's  degree  at  Harvard.  He  continues  to  be  associated  with 
the  Office  of  Naval  Research.  Roberts  has  received  his  degree  since  the  last 
Quarterly  Progress  Report,  and  has  left  to  take  a  position  at  the  Bell  Telephone 
Laboratories,  Inc. 
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Professor  M.  P.  Barnett,  who  has  been  with  us  since  1958,  has  just 
been  appointed  to  a  Readership  in  Information  Processing  at  the  Computer  Unit 
of  the  University  of  London,  a  position  in  which  he  will  be  a  senior  member  in 
one  of  the  largest  computer  installations  in  the  world.  He  will  leave  in 
September  1963  to  take  up  his  new  post.  He  will  leave  a  serious  gap  in  our 
group,  and  arrangements  have  not  yet  been  made  to  fill  his  place,  either  as 
Director  of  the  Cooperative  Computing  Laboratory  or  as  a  scientific  member 
of  the  group.  It  is  our  intention,  however,  to  continue  the  operation  of  the 
Cooperative  Computing  Laboratory  along  its  present  lines,  which  have  been  so 
soundly  originated  by  Prof.  Barnett,  and  to  continue  the  work  on  molecular  cal¬ 
culation  and  molecular  integrals.  It  is  our  hope  that  during  the  remaining 
months  of  Prof.  Barnett’s  term  here,  with  the  large  group  now  working  on  the 
programming  problem  for  the  polyatomic  molecules,  we  shall  be  able  to  bring 
this  problem  to  a  point  where  it  provides  usable  methods  for  the  quantitative 
study  of  polyatomic  molecules. 


J.  C.  Slater 
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ENERGY  BANDS  FOR  THE 


IRON  TRANSITION  SERIES 


L.  F.  Mattheiss 


I.  Introduction 

Energy  bands  have  been  calculated  for  a  majority  of  the  elements  in 

(1  2) 

the  iron  transition  series  using  the  augmented  plane  wave  (APW)  method  ’  . 

While  these  results  are  preliminary  in  nature  and  not  in  any  sense  complete, 
they  may  be  of  some  interest  to  experimentalists  and  theoreticians  who  are 
concerned  with  the  electronic  structure  of  the  transition  series  elements.  The 
present  results  represent  energy  bands  for  three  different  crystal  structures, 
with  a  variety  of  lattice  constants.  Despite  the  detailed  differences  that  are 
imposed  by  symmetry  requirements  and  variations  in  lattice  constants,  the 
results  suggest  some  interesting  and  rather  clear  cut  trends  in  the  band 
structure  of  these  elements  as  one  proceeds  through  the  transition  series. 

These  calculations  lend  some  support  to  the  rigid  band  model  for  the  transition 
series.  They  support  the  hope  that  systematic  studies  of  the  band  structure  of 
the  transition  series  elements  can  provide  useful  qualitative,  and  perhaps  quanti¬ 
tative,  information  concerning  their  electronic  structure. 

As  in  all  calculations  involving  d  electrons,  the  results  are  sensitive 
to  the  choice  of  potentials.  The  crystal  potentials  used  in  these  calculations 
were  all  constructed  in  an  analogous  manner,  and  were  approximated  by  a 
superposition  of  atomic  potentials.  The  method  involves  the  use  of  Hartree- 

(3) 

Fock  solutions  to  the  corresponding  atomic  problem  and  the  free  electron 

/4\ 

exchange  approximation'  The  details  of  this  method  for  constructing 

(5) 

approximate  crystal  potentials  have  been  described  earlier  ,  though  a  brief 
resume  is  presented  in  Section  II  of  this  paper,  along  with  other  information 
pertaining  to  this  present  series  of  calculations.  The  energy  bands  along  a 
single  line  of  symmetry  in  the  appropriate  Brillouin  zone  are  presented  in 
Section  III  for  A,  Ti,  V,  Cr,  Fe,  Co,  Ni,  Cu,  and  Zn,  while  the  last  section 
contains  a  brief  discussion  of  these  results. 
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11.  De&cription  of  the  Calculations 

In  these  calculations,  the  crystal  potential  has  been  approximated  by  a 
superposition  of  atomic  potentials.  The  coulomb  and  exchange  contributions 
to  the  crystal  potential  are  treated  separately.  An  approximate  crystal  coulomb 
potential  and  charge  density  in  a  given  atomic  cell  is  obtained  by  expanding  the 
neutral  atom  coulomb  potentials  and  charge  densities  of  neighboring  atoms  about 
the  origin,  using  Lowdin's  alpha  function  expansion^^^  keeping  only  the  i  =  0 
or  spherically  symmetric  terms  in  these  expansions.  Using  the  free  electron 
exchange  approximation,  the  exchange  potential  is  proportional  to  the  cube  root 
of  the  superimposed  atomic  charge  densities.  It  is  felt  that  by  using  this 
superimposed  charge  density,  the  principal  deficiency  in  the  free  electron  ex¬ 
change  approximation  is  compensated,  namely  that  it  falls  off  too  rapidly  at 
large  radii  in  atomic  systems. 

The  potentials  obtained  by  this  method  are  generally  rather  flat  near 
the  boundaries  of  the  atomic  cell,  at  least  in  the  case  of  metals,  so  they  are 
readily  approximated  by  a  ‘'muffin-tin"  type  potential,  as  required  by  the  APW 
method.  The  constant  value  of  the  potential  outside  the  APW  spheres  is  taken 
as  the  average  value  of  the  potential  in  this  region.  This  usually  results  in  a 
discontinuity  in  the  potential  at  the  sphere  radius  amounting  to  a  few  hundredths 
of  a  Rydberg. 

In  the  construction  of  approximate  crystal  potentials  for  transition 
series  elements,  there  is  frequently  some  ambiguity  in  choosing  the  most 
reasonable  atomic  configuration.  This  sort  of  difficulty  can  only  be  answered 
satisfactorily  by  experimental  information  and/  or  self  consistent  energy  band 
calculations.  For  the  present,  we  have  been  content  to  study  the  effect  that 
changing  the  atomic  configuration  has  on  the  band  structure.  In  addition,  there 
are  magnetic  effects  which  create  additional  complications  in  this  series  of 
elements.  For  simplicity,  all  magnetic  effects  have  been  neglected  in  these 
calculations,  and  the  crystals  have  been  assumed  to  be  non^magnetic  in 
character. 

The  lattice  constants  which  have  been  used  in  these  calculations  have 
generally  been  the  room  temperature  values  as  tabulated  in  Pearson’s  book, 

A  Handbook  of  Lattice  Spacings  and  Structures  of  Metals  and  Alloys  .  The 
exceptions  are  those  for  A  and  Zn.  The  lattice  constant  for  A  is  the  low  temp¬ 
erature  value  obtained  by  Dobbs  and  Jones In  the  case  of  Zn,  Harrison^*^^ 
has  extrapolated  the  room  temperature  lattice  constants  to  low  temperatures 
since  the  results  are  expected  to  be  sensitive  to  the  choice  of  c/  a  ratio.  For 
purposes  of  comparison,  his  values  have  been  used  in  these  calculations. 

Table  1  contains  a  summary  of  the  elements  considered  in  these  calculations, 
their  structures,  the  values  of  the  lattice  constants,  and  the  assumed  atomic 
configurations. 
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Table  1 

In  this  table,  we  list  the  elements  of  the  iron  transition  series,  their  structures, 
the  lattice  constants  used  in  this  series  of  calculations  (in  atomic  units),  and  the 
assumed  atomic  configurations. 


Element 

Structure 

a  (au) 

c  (au) 

Configuration 

A 

fee 

10. 0346 

(3s)^(3p)^ 

K 

Ca 

Sc 

Ti 

hep 

5.  5755 

8.8503 

(3d)^4s)^ 

V 

bcc 

5.  7225 

(3d)^(4s)^ 

Cr 

bcc 

5.  4512 

(3d)^4s)^ 

Mn 

Fe 

bcc 

5. 4168 

(3d)'^(4s)^ 

Co 

fee 

6. 6975 

(3d)®(4s)^ 

Ni 

fee 

6. 6590 

(3d)^(4s)^ 

Cu 

fee 

6.  8309 

(3d)^°(4s)^ 

Zn 

hep 

5. 0120 

9. 1453 

(3d)^°(4s)^ 

III.  Results 

The  principal  results  of  these  calculations  are  presented  in  Figure  1. 
These  results  represent  plots  of  energy  as  a  function  of  wave  vector  along  lines 
of  symmetry  from  the  center  to  a  boundary  of  the  appropriate  Brillouin  zones. 
For  the  face  centered  cubic  structures  (A,  Co,  Ni,  and  Cu),  the  bands  are 
plotted  from  F  along  the  A  direction  to  the  point  X,  using  the  notation  of 
Bouckaert,  Smoluchowski,  and  Wigner^^^^  In  the  body  centered  cubic  structure 
(V,  Cr,  and  Fe),  the  bands  are  plotted  from  T  along  the  A  direction  to  the  point 
H.  Finally,  in  the  hexagonal  close  packed  structure  (Ti  and  Zn),  they  start  at 

r  and  proceed  along  the  line  T  in  the  k  =0  plane  which  terminates  at  the  point 

^  (11) 

K,  one  of  the  vertices  of  the  hexagon  (in  the  notation  of  Herring' 

The  energy  is  in  Rydbergs  and  the  wave  vectors  for  the  different 
elements  are  drawn  to  scale  for  purposes  of  comparison.  The  horizontal 
dashed  lines  represent  rough  estimates  of  the  Fermi  energy  for  each  element. 
For  simplicity,  some  of  the  more  highly  excited  bands  have  been  omitted  in 
some  cases,  particularly  in  the  face  centered  cubic  structure,  or  in  other 
situations,  they  have  been  sketched  in  by  dashed  lines. 
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Figure  1 


Energy  bands  for  A,  Ti,  V,  Cr,  Fe,  Co,  Ni,  Cu,  and  Zn  as  a  function  of 
wave  vector  along  a  line  of  symmetry  in  the  appropriate  Brillouin  zones. 
For  the  face  centered  cubic  structure,  the  bands  are  plotted  from  F  along 
^  to  the  point  X.  For  the  body  centered  cubic  structure,  they  are  plotted 


ENERGY  BANDS  FOR  THE  IRON  TFlANSITION  SERIES 


9. 


1.0 


E 

(Rydbergs) 


0.5 


0 


Figure  1  (cont. ) 

from  r  along  A  to  the  point  H.  Finally,  for  the  hexagonal  close  packed 
structure,  they  are  plotted  from  V  along  T  to  the  point  K.  The  energies 
are  in  Rydbergs  and  the  wave  vectors  are  in  atomic  units. 
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In  Figure  2,  it  is  shown  what  effect  varying  the  atomic  configuration 
has  on  the  band  structure  of  a  typical  element,  namely  vanadium.  The  bands 
to  the  right  are  the  ones  shown  in  Figure  1  for  vanadium;  the  ones  to  the  left 
are  those  obtained  from  a  potential  which  results  from  an  atomic  configuration 
containing  an  additional  4s  electron  and  one  less  3d  electron.  In  general,  this 
results  in  a  narrowing  of  the  3d  band  and  a  decrease  in  the  energy  separation 
between  the  top  of  the  3d  band  and  the  bottom  of  the  4s-4p  bands. 

IV.  Discussion 

In  the  simplified  picture  of  the  energy  bands  for  the  iron  transition 
series  elements,  one  finds  a  narrow  3d  band  in  the  midst  of  a  rather  broad 
4s-4p  band.  The  width  of  the  3d  band  and  especially  its  position  relative  to 
the  bottom  of  the  4s”4p  band  depend  rather  critically  on  the  potential.  Never¬ 
theless,  the  results  of  Figure  1  exhibit  a  reasonably  smooth  variation  from 
element  to  element,  especially  for  those  substances  having  the  same  crystal 
structure.  This  seems  to  lend  some  support  to  the  rigid  band  model  for  the 
transition  series  elements,  an  approximation  which  has  been  of  considerable 
value  in  understanding  the  electronic  properties  of  these  elements  and  their 
alloys. 

There  is  a  gradual  narrowing  of  the  3d  band  as  one  progresses  through 

the  series.  This  effect  was  discussed  by  Slater  in  order  to  explain  the 

(12) 

occurrence  of  ferromagnetism  in  the  latter  part  of  the  series'  ’ .  In  going 
from  Cu  to  Zn,  the  3d  band  suddenly  drops  about  0.  5  Rydberg  below  the  bottom 
of  the  4s -4p  bands,  and  its  width  decreases  to  less  than  0.  1  Rydberg.  As  a 
result,  the  energy  bands  for  zinc  are  very  free  electron-like.  For  those  ele¬ 
ments  where  the  3d  band  falls  in  the  middle  of  the  4s-4p  bands,  the  inter¬ 
actions  between  states  having  the  same  symmetry  causes  considerable  modi¬ 
fication  to  the  free  electron  bands,  though  at  points  of  symmetry,  the  effect 
is  sometimes  small.  The  bands  for  Ti  and  Zn  demonstrate  this  effect  nicely. 

The  results  of  Figure  2  emphasize  the  uncertainty  which  is  inherent 
in  any  energy  band  calculation  for  a  transition  series  element.  These  uncer¬ 
tainties  have  been  pointed  out  previously  in  the  literature,  particularly  by 
(13) 

Callaway'  These  difficulties  can  only  be  cleared  up  satisfactorily  with  the 
aid  of  more  detailed  experimental  information  regarding  the  band  structure  of 
these  elements  in  addition  to  self  consistent  energy  band  calculations. 

The  results  presented  here  are  not  complete  enough  to  permit  de¬ 
tailed  comparisons  to  be  made  with  experiment  or  a  discussion  of  the  resulting 
Fermi  surfaces.  However,  there  are  some  striking  similarities  between  the 
energy  bands  shown  in  Figure  1  and  the  results  obtained  by  earlier  calculations. 

In  particular,  there  is  good  qualitative  agreement  between  the  Cu  results 

(14)  (15) 

shown  on  Figure  1  and  the  bands  calculated  by  Segall  ’  and  also  by  Burdick 
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Energy  bands  for  V  using  potentials  obtained  from  two  different  atomic 

configurations.  The  bands  to  the  left  resulted  from  a  configuration  of 

3  2  4  1 

(3d)  (4d)  ,  while  those  to  the  right  involved  a  (3d)  (4s)  configuration. 
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Similarly,  the  results  for  Fe  are  in  good  agreement  with  the  published  results 
of  Wood^^\  The  agreement  for  A  with  the  results  of  calculations  by  Knox  and 
Bassani^^^^  is  good,  and  has  been  described  previously^^\  though  there  are 
differences  in  the  ordering  of  levels.  In  the  case  of  Cr,  it  is  difficult  to  com¬ 
pare  the  present  results  with  those  of  earlier  calculations  by  Asdente  and 
Friedel^^^^  since  they  neglect  the  interactions  between  the  3d  band  and  the 
4s -4p  bands.  However,  the  present  results  do  justify,  to  some  extent,  the 
treatment  of  Cr  by  Lomer^^^^,  who  used  the  results  of  Wood's  iron  calculations 
to  discuss  the  energy  bands  in  antiferromagnetic  Cr. 

In  the  case  of  Zn,  the  ordering  of  levels  is  identical  with  that  obtained 
(g\ 

by  Harrison'  \  This  ordering  differs  from  that  obtained  in  earlier  calculations 

(19) 

for  hexagonal  close  packed  metals  by  Herring  and  Hill  for  Be'  '  and  Falicov 

for  This  change  in  ordering  might  be  due  to  the  presence  of  an 

occupied  3d  band  just  below  the  4s -4p  bands.  Finally,  the  results  for  Ti  agree 

(21 ) 

qualitatively  with  those  obtained  by  Altmann  and  Bradley' 
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ENERGY  BANDS  IN  GALLIUM 


J.  H.  Wood 


The  figures  shown  in  the  last  progress  report  should  be  disregarded 
as  they  are  considerably  in  error.  We  now  believe  that  the  structure  of  the 
5-6  Fermi  surface  (which  appears  to  be  the  important  one  experimentally) 
resembles  the  corresponding  free  electron  surface  shown  in  figure  i4{a)  of 
the  paper  by  Reed  and  Marcus The  main  modifications  (which  are  going 
to  be  considerable)  to  that  surface  are  in  the  regions  of  the  reflection  plane 
perpendecular  to  the  axis  (in  that  figure)  --  this  is  the  region  in  which  the 
free  electron  surface  allows  an  "exit"  from  the  zone,  in  the  vicinity  of  the 
vertices  of  the  hexagon.  We  are  currently  concentrating  our  calculations  in 
that  region  of  k-space  in  order  to  ascertain  if  an  exit  is  still  possible  in  this 
region. 

Some  work  has  also  been  going  into  freeing  the  APW  programs  from 
the  muffin  tin  type  potential.  This  type  of  potential  is  somewhat  restrictive 
whenever  the  volume  of  space  outside  the  muffins  becomes  a  sizeable  fraction 
of  the  total  cell  volume.  We  do  not  intend  to  immediately  change  over  tO'  this 
for  the  case  of  gallium,  even  though  some  60  percent  of  the  cell  volume  is 
outside  the  APW  spheres  and  is  being  represented  by  a  constant  value  for  the 
potential.  It  seems  worthwhile  to  push  the  present  calculations  to  a  conclusion 
to  see  if  they  will  give  us  something  sensible  in  the  way  of  a  Fermi  surface 
W'hich  would  be  of  use  to  the  experimental  people. 

We  also  hope  to  soon  begin  testing  the  program  for  the  diamond  type 
structure . 


Reference 


1.  W.  A.  Reed  and  J.  A.  Marcus,  Phys.  Rev.  126,  1298  (1962) 
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The  space  group  having  the  point  group  and 

the  primitive  tetragonal  Bravais  lattice,  is  found  in  the  rutile  structure, 
occurring  in  TiO^  and  similar  compounds.  The  three  fundamental  vectors 
tf,  of  the  primitive  unit  cell  are  given  by  t^  =  ai,  t^  =  aj,  ^  =  c  k, 

and  the  vectors  b^,  b^,  of  the  reciprocal  lattice  are  bj  =  i/  a,  b^  =  j/  a, 
b^  =  k/  c.  The  Brillouin  zone,  and  the  symmetry  points  in  it,  are  given  in 
Fig.  1.  The  positions  of  the  atoms  in  the  rutile  structure  are  as  follows: 

Titanium  at  the  origin  of  the  cell,  and  at  the  point  for  which  x/  a,  y/  a,  z/  c 
are  I/2.  l/2,  l/2. 

Oxygen  at  =t(uuO;  u+l/ 2,  1/  2-u,  l/  2),.  where  for  TiO^  the  parameter  u  is 
0.  31. 

The  16  operations  of  the  point  group  are  Xq.  x^.  y^. 

Y^.  X^,  x;^,  X^.  Y;,,  Y;^j,  Y^,  defined  by  Ers.  (6)  .  .  .  (9). 
Reference  1.  That  is,  we  have 

X  ij^(r,  <j>,  z)  =  <|)  +  -5^,  z) 

^  N 

y  ij>(r,  <|),  z)  =  t(r,  -<|>  +  z) 

^  N 

X* ‘^(r,  c|>,  z)  =  i^r,  (J>  +  ,  -z) 

^  N 

Y' ^{r,  <j),  z)  =  +(r,  -<j>  +  -z) 

^  M 


(1) 
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where  N  =  4.  For  the  present  purpose  it  is  more  convenient  to  express  the 
operations  in  rectangular  coordinates.  We  then  have 


Xq+Cx,  y,  z)  =  i|/(x,  y,  z) 
Xj+(x,  y,z)  =  4'(-y.x,  z) 
X_j+(x,  y,  z)  =  i|j(y,  -X,  z) 
X2+(x,  y,  z)  =  i>(-x,  -y,  z) 
Yq+(x.  y,  z)  =  i}/(x,  -y,  z) 
Yji|j(x,y,z)  =  i|/(y,  X,  z) 

Y_  j>|j(x,  y,  z)  =  4/(-y,  -x,  z) 
Y,>j'(x,  y,  z)  =  iW-x,  y,  z) 


Xq^Kx,  y,  z)  =  i|j(x,y,  -z) 
X^v^x,  y,  z)  =  ^(-y.x.  -z) 
X;_  j^|j(x,  y,  z)  =  4/(y,  -X,  -z) 
X^i|j(x,  y,  z)  =  v)j(-x,  -y,  -z) 
Y^>Jj(x,  y,  z)  =  4/(x,  -y,  -z) 
Y^Kx,  y,  z)  =  i|j(y,  x,  -z) 
Yl^jWx.y,  z)  =  i|j{-y,  -x,  -z) 
Y''l'(x,  y,  z)  =  >l'(-x,y.  -z) 


(2) 

The  multiplication  table  for  the  point  group  is  as  given  in  Eq.  (5),  Ref.  1, 
namely 


,  X„Y  =  Y  .  Y  X  =  Y  ,  Y^Y  =  X  ... 

q  p  q+p  q  p  -q+p  q  p  p+q’  q  p  -q+p  (3) 


In  this  expression,  if  the  subscript  of  the  product  function  does  not  have  one 

of  the  allowed  values  0,  ±1,  2,  we  are  to  add  or  subtract  integral  multiples 

14 

of  4  to  secure  a  value  within  this  range.  In  the  space  group  we  have  no 

non-primitive  translation  associated  with  operators  with  even  subscripts,  but 
the  non-primitive  translation  ^  associated  with  all  operators  with 

odd  subscripts. 

The  effect  of  the  operators  of  Eq.  (2)  on  a  plane  wave 
exp  Zui  [(h^+p^)4+(h2+p2)'n+(h2+P2)^3]  »  where  the  vector  position  7  is  given 
^  given  by  Eq.  (4)  below: 

Xq1|j  =  exp  Z-ni  [  (hj+Pj)g+(h2+P2)-n+(h3+P2)C  ] 

Xj+  =  exp  2i7i  [  (h2+P2)5-(hj+Pj)Ti+(h2+P2)C  ] 

X_ji<J  =  exp  Z-ni  [  -(h2+P2)S+(hj+Pj)T|+(h2+P2)4] 

X2'l'  »  exp  2iri  [  -(hj+Pj)^-(h2+P2)i1+(h3+P3)C  ] 

Yq+  =  exp  2iii  [  (hj+Pj)g-(h2+P2)r|+(h3+P3)^  ] 

Yjij/  =  exp  2iri  [  (h2+P2)l+(hj+Pj)T|+(h2+P2)^  ] 

Y_ji|j  =  exp  2iri  [ -(h2+P2)^-(hj+Pj)T|+(h2+P3)4  ] 

Y2't'  =  exp  Z-ni  [  -(hj+Pj)^+(h2+P2)h+(h3+P3)&  ] 


(4) 
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For  the  primed  operations,  the  sign  of  h^+p^  is  to  be  changed.  An  additional 
factor  exp  'rri(h^+p^+h24-p2+h^+p^)  occurs  for  the  operations  with  odd  sub¬ 
scripts  in  the  space  group,  on  account  of  the  non-primitive  translations.  In 
the  multiplication  table  of  the  space  group,  there  is  an  additional  factor, 
arising  from  the  non-primitive  translation,  if  the  first  operation  in  the  pro¬ 
duct  has  an  odd  subscript.  There  is  no  such  factor  if  the  first  operation  has 
an  even  subscript.  The  factor  is  independent  of  which  odd-subscript  operator 
we  have,  and  is  as  follows,  where  we  first  enumerate  the  second  operation, 
then  the  factor: 


: 

1 

0 

Xi  : 

^2^i(P2+P3) 

X  ,  : 

-1 

X3  : 

e'2’ri(Pi+P2) 

Yq  = 

= 

g2TTi(Pj+p2+p. 

Y-r 

e2^iP3 

^2  = 

e-2iTiPi 

e-2^iP3 

e2TviPj 

g-2Tri(Pj+P2+P3) 

g-2ni(p2+P3) 

e2Tri(Pi+P2) 


1 

g-2ui(pj+P3) 


(5) 


We  now  give  in  Tables  1  ...  14  the  matrix  elements  and 
characters  for  the  various  symmetry  directions,  and  in.  Table  15  the  com¬ 
patibility  relations. 
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Table  1 

Matrix  elements  and  characters  at  point  F,  =  0.  Only  the  matrix 

elements  for  the  unprimed  operators  are  tabulated.  For  the  primed  operators, 
the  elements  are  the  same  as  for  the  unprimed  for  the  representations  with 
superscript  +,  and  the  negatives  of  these  values  for  the  representations  with 
superscript 


Table  2 


Matrix  elements  at  points  A,  p,  =  pi  p-  =  p«  =  0 
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Table  3 

Matrix  elements  at  point  X,  boundary  of  zone  along  direction  A.  Here 

Pi  =  1/^'  Pz  "  P3  "  °- 


^0 

^2 

^0 

^2 

^0 

^2 

Y' 

0 

^2 

^1 

1 

1 

1 

1 

1 

1 

1 

1 

^2 

1 

-1 

1 

-1 

1 

-1 

1 

-1 

X3 

1 

-1 

-1 

1 

1 

-1 

1 

^4 

1 

1 

1 

1 

-1 

-1 

-1 

-1 

^5 

1 

-1 

-1 

1 

-1 

1 

1 

-1 

^6 

1 

-1 

1 

-1 

-1 

1 

-1 

1 

X7 

1 

1 

-1 

-1 

1 

1 

-1 

00 

1 

1 

-1 

-1 

-1 

-1 

1 

1 

Table  4 

Matrix  elements  and  characters  for  points  A,  p,  =  Po  =  P- 


2 


0 


0 


0 


0 
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Table  6 

Matrix  elements  for  points  S,  =  p,  =  0 


^0 

^1 

^'0 

Y' 

1 

1 

eZuip 

1 

eZTtip 

^2 

1 

.e2^iP 

-1 

e2uip 

^3 

1 

eZuip 

-1 

_e2,rip 

^4 

1 

_e2^ip 

1 

_e2uip 

Table  7 

Matrix  elements  and  characters  for  point  M,  boundary  of  zone  along  direction 
2.  Here  Pj  P2  ~  P3  “  table  below  the  matrix  elements  for 

the  primed  operations  are  equal  to  those  for  the  unprimed  for  superscript  +, 
and  the  negative  of  these  values  for  superscript  ^  . 


^0 

^-1 

^2 

^0 

Yi 

Y-1 

^2 

1 

1 

-1 

-1 

1 

1 

-1 

-1 

1 

1 

-1 

-1 

-1 

-1 

1 

1 

1 

-1 

1 

-1 

1 

-1 

1 

-1 

M  ± 

4 

1 

-1 

1 

-1 

-1 

1 

-1 

1 

(^5=^)1! 

1 

0 

0 

1 

0 

1 

1 

0 

(^5=^)2! 

0 

1 

1 

0 

1 

0 

0 

1 

(M5*)i2 

0 

-1 

-1 

0 

1 

0 

0 

1 

<M5*)22 

1 

0 

0 

1 

0 

-1 

”1 

0 

X  (Mgi) 

2 

0 

0 

2 

0 

0 

0 

0 
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Table  8 

Matrix  elements  for  directions  W,  =  1/  2,  =  0,  =  p.  For  point  X, 

p  =  0,  For  point  R,  p  =  l/  2. 


^0 

^2 

^0 

^2 

Wi 

1 

1 

1 

1 

^2 

1 

1 

-1 

-1 

W3 

1 

-1 

1 

-1 

W4 

1 

-1 

-1 

1 

Table  9 

Matrix  elements  for  direction  R,  =  P3  “  P2  ”  ^ 


^0 

^2 

^0 

^2 

^0 

^2 

1 

1 

1 

1 

1 

1 

1 

1 

«2 

1 

-1 

1 

-1 

1 

-1 

1 

-1 

R3 

1 

-1 

-1 

1 

•1 

-1 

-1 

1 

R4 

1 

1 

1 

1 

-1 

-1 

-1 

-1 

R5 

1 

-1 

1 

-1 

1 

1 

-1 

1 

-1 

1 

-1 

-1 

1 

-1 

1 

R? 

1 

1 

-1 

-1 

1 

1 

-1 

-1 

Rr 

1 

1 

-1 

-1 

-1 

-1 

1 

1 
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Table  10 

Matrix  elements  and  characters  at  points  V,  joining  M  and  A.  Here 


Pj  =P2  =  1/2. 

P3  =P. 

At  M, 

*0 

il 

0 

at  A,  p 

=  1/2. 

x„ 

X, 

X  , 

x^ 

Y, 

Y  , 

Y^ 

0 

1 

-1 

2 

0 

1 

-1 

2 

Vi 

1 

e^iP 

irip 

-e 

-1 

1 

e^iP 

-1 

^2 

1 

irip 

e 

_e^iP 

-1 

-1 

-rrip 

-e 

irip 
e  ^ 

1 

V3 

1 

^irip 

-1 

1 

e^iP 

-1 

^4 

1 

_e^iP 

e^iP 

-1 

-1 

e^iP 

.e^iP 

1 

<Vll 

1 

0 

0 

1 

1 

0 

0 

1 

^^5^21 

0 

gtrip 

e’fiP 

0 

0 

_e’fiP 

.e’fiP 

0 

<V5>12 

0 

.e’fiP 

0 

0 

.e^ip 

0 

<^5)22 

1 

0 

0 

1 

-1 

0 

0 

-1 

x(V 

2 

0 

0 

2 

0 

0 

0 

0 
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SPACE  GROUP  D^^(P42/  mnm) 


Table  11 


Matrix  elements  and  characters  for  point  A, 


Pj  =  P2  =  P3  =  1/  2 


^0 

^1 

x-1 

^2 

^0 

Yi 

Y-1 

^2 

^0 

^2 

Y’ 

Y’l 

YLi 

(Aj)  jj 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

*^1*  21 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

0 

0 

0 

0 

<^1)12 

0 

0 

0 

0 

1 

-1 

-1 

1 

1 

-1 

-1 

1 

0 

0 

0 

0 

*■^1^22 

1 

-1 

-1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

-1 

-1 

1 

XtAj) 

2 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

2 

*■^2^1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

-1 

-1 

-1 

<^2)21 

0 

0 

0 

0 

1 

1 

1 

1 

-1 

-1 

-1 

-1 

0 

0 

0 

0 

*■^2^2 

0 

0 

0 

0 

1 

-1 

-1 

1 

-1 

1 

1 

-1 

0 

0 

0 

0 

*■^2^22 

1 

-1 

-1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

1 

1 

-1 

XIA^) 

2 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

-2 

0 

0 

-2 

*■^3*11 

1 

I 

0 

0 

-1 

1 

0 

0 

-1 

1 

0 

0 

-1 

1 

0 

0 

-1 

<^>21 

0 

1 

-1 

0 

0 

1 

-1 

0 

0 

-1 

1 

0 

0 

-1 

1 

0 

*^3^12 

0 

-1 

1 

0 

0 

-1 

1 

0 

0 

-1 

1 

0 

0 

-1 

1 

0 

<^3)22 

1 

0 

0 

-1 

1 

0 

0 

-1 

-1 

0 

0 

1 

-1 

0 

0 

1 

x(A3) 

2 

0 

0 

-2 

2 

0 

0 

-2 

0 

0 

0 

0 

0 

0 

0 

0 

(•^4^1 

1 

0 

0 

-1 

-1 

0 

0 

1 

-1 

0 

0 

1 

1 

0 

0 

-1 

*■^4^21 

0 

1 

0 

0 

-1 

1 

0 

0 

1 

-1 

0 

0 

-1 

1 

0 

*■^4^2 

0 

-1 

1 

0 

0 

1 

-1 

0 

0 

1 

’1 

0 

0 

-1 

1 

0 

<^4)22 

1 

0 

0 

-1 

-1 

0 

0 

1 

1 

0 

0 

-1 

-1- 

0 

0 

1 

XIA^) 

2 

0 

0 

-2 

-2 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 
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Table  12 

Matrix  elemeats  at  points  U,  connecting  Z  and  R.  Here  =  p,  ~  0, 
p^  =  1/2.  At  Z,  p  =  0;  at  R,  p  =  i/2. 


^0 

Yo 

^'0 

n 

1 

1 

1 

1 

U2 

1 

-1 

-1 

1 

U3 

1 

-1 

1 

-1 

U4 

1 

1 

-1 

~i 

Matrix  elements  at  points  T, 

Table  13 

connecting  R  and  A. 

Here  p^ 

p^  =  p.  At  R,  p  =  0; 

at  A,  p 

X 

0 

=  1/2. 

^2 

1 

1 

1 

1 

Tz 

1 

-1 

-1 

1 

T3 

1 

-1 

1 

-1 

T. 

1 

1 

-1 

-1 

Table  14 

Matrix  elements  and  characters  at  points  S,  connecting  Z  and  A,  Here 
Pj  =  P2  “  P*  P3  ~  At  Z,  p  =  0;  at  A,  p  =  1/  2. 
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SPACE  GROUP  D^^(P42/mnm) 


Table  15 

Compatibility  relations 


1 

^2 

^2 

p  + 

3 

^3 

^  4 

^4  ^5  ^5 

^4 

^3 

^2 

^1 

^4 

^3 

^2  ^1^  ^2^4 

^2  ' 

■^2 

A3 

A3 

^4 

^>1  Ac  Ac 

4  5  5 

^3 

^4 

^2 

^4 

^2 

^1 

Z3  Z2Z3 

^4 

^5 

^6 

X7 

^8 

^1 

^3 

^4 

^2 

^4 

^3 

^2 

Y3 

^4 

^4 

^2 

Y3 

^2 

Wi 

^^3 

W4 

Wj 

W4 

W3 

^2 

^2 

^2 

^3 

^4 

A1A3 

A^A 

4 

^5 

^5 

U1U4 

^2^ 

3 

U1U2 

U3U4 

Si 

Si 

Si 

m; 

M- 

^2 

^2 

M3 

M3 

M+ 

m'  m^  m" 

4  5  5 

^4 

^2 

^3 

^1 

^3 

^4 

^2  ^1^4  ^2^: 

^3 

^2 

Yl 

Y4 

Y3 

^2 

Yl 

Y4  Y1Y3  Y3Y, 

Vl 

^2 

^2 

V3 

^3 

^4 

^4  ^5  ^5 

«1 

«2 

R3 

«4 

^5 

«6 

«7 

00 

^1 

Ul 

U3 

U4 

^2 

U4 

U3 

U2 

Wi 

W, 

3 

W4 

Wi 

W4 

W3 

^2 

^3 

^1 

^4 

T4 

^2 

'^3 

^2 
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SPACE  GROUP  D^^(P42/mnm) 
Table  15  (continued) 


^2 

A3 

^4 

Si 

Si 

Si 

^1^2 

^3^4 

T2T3 

T1T4 

^5 

^5 

^1^3 

^^2^4 
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SYMMETRY  PROPERTIES  IN 
C3.(R3),  D^(R32),  D^^(RJm) 


THE  SPACE  GROUPS 
,  AND  Dj^iRlc) 


John,  C.  Slater 


The  space  groups  considered  here  have  the  rhombohedral  or  trigonal 

Bravais  lattice,  and  can  be  conveniently  considered  together.  A  good  many 

important  elements  and  compounds  crystallize  in  one  or  another  of  these 

space  groups.  Thus  the  ilmenite  structure  has  the  space  group  C?.{R3),  the 
7 

AIF^  structure  has  D2(R32),  the  arsenic  structure,  the  CdCl2  structure,  ana 
the  LaOF  structure  have  the  space  group  D|^(R3m)  ,  and  the  calcite  and 
corundum  structures  have  the  space  group  D°^(R3c).  Positions  of  atoms  in 
the  structures  listed  above  are  given  in  Table  1. 

The  three  fundamental  vectors  Tj,  of  the  primitive  unit  cell  are 

of  equal  magnitude,  and  make  equal  angles  of  a  with  each  other.  In  rectangu¬ 
lar  coordinates  they  are  expressed  as 

Tj  =  sr+  r1c,  t*^  =  +  ric,  t^  =  (1) 


^  — 2 — ' — 2” 

The  length  of  the  vectors  t^,  T^,  is  s/s^+r  ,  and  the  angle  a  is  given  by  the 
equation 


cos  a 


-s 


—Z — z 


s  +r 


2 

r 


(2) 


from  which  we  can  find  s  and  r  when  the  length  of  the  vector  and  the  angle 
a  are  given,  as  is  usual  in  crystallographic  work.  The  vectors  b^, 
of  the  reciprocal  lattice  are  given  by 


^  =  s'r+  r’?,  =  —  (-T  +  +  r'lT,  ?  =  —  ‘  +  r'iT  (3) 

^  ^  1 
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SPACE  GROUPS  Cj.,  d],  AND 

where  s'  =  2/ 3s,  r'  =  l/3r  (4) 

From  Eq.  (3)  we  see  that  the  reciprocal  lattice  is  also  rhombohedral.  The 
Brillouin  zone  is  shown  in  Fig.  1,  with  convenient  notations  (partly  following 
Koster^)  for  the  symmetry  directions.  The  inclined  hexagonal  faces  are  per¬ 
pendicular  bisectors  of  vectors  rectangular  faces  are 

perpendicular  bisectors  of  the  vectors  ±  ^  and  the 

top  and  bottom  faces  are  perpenticular  bisectors  of  the  vectors  ± 

The  height  of  the  Brillouin  zone  is  3r',  and  the  dividing  lines  between  the 
rectangular  faces  and  the  inclined  hexagonal  faces  come  at  distances  r'  and 
2r'  from  the  top  and  bottom  faces.  The  symmetry  points  indicated  by  A  and 
B  in  Fig.  1  are  the  intersections  of  the  vectors  **5^,  etc.  ,  and  ±  etc.  , 

with  the  faces  of  the  Brillouin  zone. 

The  point  groups  and  D^,  which  are  found  in  the  first  two  space 
groups  considered,  may  be  taken  as  special  cases  of  with  the  re¬ 
maining  two.  Consequently  we  shall  carry  through  the  case  of  and  shall 

specialize  for  the  other  cases.  The  operations  of  the  point  group  in 

cylindrical  coordinates,  may  be  taken  to  be  defined  by 

XgiKr,  4),  z)  =+(r,  4>,  z) 

4),  z)  =  4j(r,  4>ij,  -z) 

=  4'(r.  4>=t^.z) 

(r,  <j),  z)  =  i|i(r,  (|>+1T,  -z) 

Y^iKr,  <}>,  z)  =+(r, -<|>,  z) 

Y^j  +  (r,  4),  z)  =  -4)±j,  -z) 

(r,  4),  z)  =  4>(r,  -4)+Tr,  -z)  (5) 

The  operations  of  may  be  described  in  terms  of  Eq.  (5)  as  the  X  opera¬ 
tions  listed  there;  the  operations  of  may  be  described  as  Xq,  X^2» 

Y^.  In  Table  2  we  describe  the  operations  of  in  terms  of  their  effect  on 
a  function  of  the  vector  7  expressed  in  the  form 

7  =  +  tiT^  +  (6) 

In  Table  3  we  show  the  effect  of  these  operations  on  a  plane  wave  of  the  form 
=  exp  2iri[  (hj+p^)S^^  +  (h2+P2)S'2+(h2+p^)S2  ] '  C^) 

where  the  reduced  wave  vector  is  2Tr{p^^^+p^^+p^^),  and  the  h's  are  integers. 
These  operations  combine  according  to  the  multiplication  table 


31. 


Fig.  I  Brillouin  zone  for  the  rhombohedral  Bravais 
lattice,  with  notation  for  symmetry  directions. 
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SPACE  GROUPS  C3.,  Dy  AND  D 


'3d 


X  X 

=  X  . 

q  p 

q+p 

X  Y 

=  Y  . 

q  p 

-q+p 

Y  X 

=  Y  . 

p  q 

p+q 

Y  Y 

=  X  ^ 

q  p 

-q+p 

in  which  if  any  subscript  lies  outside  the  range  from  -2  to  3,  we  are  to  add  or 
subtract  integral  multiples  of  6  to  bring  it  within  this  range. 

The  space  groups  under  consideration  are  symmorphic,  except  for 
the  case  of  D^^^RSc).  Consequently  the  multiplication  table  as  just  stated  can 
be  used  for  the  space  groups  as  well  as  for  the  point  groups.  For  the  case  of 
D^^(R3c),  a  non-primitive  translation  (l/^2)(t^+T2+^3)  is  associated  with  each 
of  the  Y  operations;  no  non-primitive  translation  is  associated  with  the  X's. 
Consequently  the  effect  of  one  of  the  Y  operations  on  a  plane  wave  is  found 
from  Table  3  by  multiplying  the  value  given  by  that  table  by  the  factor 


^iTi(hj+p^+h2+P2+h3+p2) 


(8) 


The  effect  of  these  non-primitive  translations  on  the  multiplication  table  for 
the  space  group  {R3c)  is  found  by  using  Eq.  (7),  and  multiplying  the  result 
by  an  extra  factor  as  follows: 

Any  Y  operating  on  an  X  with  an  odd  subscript, 
factor 

Any  Y  operating  on  a  Y  with  an  even  subscript, 

factor  e^"^<Pl^P2^P3>  (9) 

In  all  other  cases  there  is  no  extra  factor. 

We  now  give  in  the  remaining* tables  the  matrix  elements  for  the 
various  irreducible  representations  of  the  four  space  groups  under  considera¬ 
tion,  and  the  compatibility  relations.  Since  there  is  no  generally  accepted 
convention  for  the  naming  of  the  irreducible  representations  of  these  space 
groups,  we  have  introduced  names  for  them. 
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Table  1 

Positions  of  Atoms  in  Structures  Having  Space  Groups  Under  Discussion 


Ilmenite  Structure  (R^).  FeTiO^*  Iron  atoms  at  positions  5  -  u,  ti  =  u, 

C  =  u  and  5  =  -u,  t)  =  -u,  t  =  -u,  where  u  =  0,358;  titanium  atoms  for  ^  =  v, 
r|  =  V,  ^  =  V,  and  i  =  -v,  t)  =  “V,  ^  =  -v,  where  v  =  0.  142;  and  oxygens  at 
general  positions  where  i  =  0.555,  t)  =  -0.055,  and  ^  = 

0.250.  Other  similar  compounds  have  similar  parameters. 


Y  _ 

AIF^  Structure  D^(R32).  Aluminum  atoms  at  uuu  and  uuu  (that  is,  |  =  u, 
=  u,  C  =  u  and  ^  =  -u,  r|  =  -u,  i  =  -u)  where  u  =  0.  237.  One  type  of 
fluorine  at  Ovv,  vOv,  vvO,  where  v  =  0.430.  Another  type  of  fluorine  at 
l/2,  ww;  w,  l/2,  w;  ww,  l/ 2,  where  w  =  0.  070.  Other  similar  compounds 
have  similar  parameters. 


5  _ 

Arsenic  Structure  (RTm).  Arsenic  atoms  at  uuu  and  uuu,  where 

u  =  0.226  for  arsenic. 


CdCl^  Structure  (RTm).  Cadmium  atom  at  origin,  chlorine  atoms  at 

uuu  and  uuu,  where  u  is  approximately  l/4. 

LaOF  Structure  (RTm).  Lanthanum  atoms  at  uuu,  uuu,  where  u  =  0.  242. 

Fluorines  in  similar  positions  with  u  =  0.  122;  oxygens  in  similar  positions, 
with  u  «  0.370. 

Calcite  Structure  (RTc)  CaCO^.  Calcium  atoms  at  origin,  and 

1/2,  1/2,  1/2.  Carbons  at  ^l/ 4,  l/ 4,  l/ 4).  Oxygens  at  ±(l/ 4-u,  l/ 4+u,  l/  4); 

1/4+u,  1/ 4,  1/ 4-u;  1/ 4,  1/ 4-u,  l/  4-fu),  where  u  =  0.  243. 

Corundum  Structure  (RTc)  Al^O^-  Aluminum  atoms  at 

i(l/4+u,  1/4+u,  l/4-fu;  1/4-u,  1/4-u,  1/4-u),  with  u  =  0.  105.  Oxygens  at 

±(1/4-v,  1/4+v,  1/4;  l/4'fv,  1/4,  1/4-v;  l/ 4,  l/ 4-v,  l/ 4+v),  with  v  =  0.303, 
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Table  2 

Operations  of  the  Point  Group  D^^ 

XqI|j(7)  = 

=  i'(-ut;-cr2-er3) 

x_jij;(r)  =  'i-(-5rj-er2-ur3) 

X2>i'(r5  =  ij^icrj+gr^+Tirj) 

X_2‘t'(r)  =  4'(Tirj+4t2+^  t^) 

Xj+ir)  =  ^(-erj-Tir^-sr^) 

=  >i'(irjfCt^+Tit^) 

Yj+{7)  = +(-ctj-Tit^-et^) 

y_jiP{7)  =  +(-nrj-et^-^r3) 

Y.2‘f'(r)  =  <W5rj+Tit^+5t^) 

Y3+(r)  =  ^'l-eTj-Cr^-nt^) 


Table  3 

Effect  of  Operators  of  on  a  Plane  Wave 

XqI|/  =  exp  2Tri  [(hj+Pj)54{h2+p2)T)+(h3+p3)& J 
Xjil*  =  exp  Eiri  [-(h3+P3)^-(hj+pj)Ti~(h2+P2)t  ] 
X_j+  =  exp  2^i  [-(h2+p2)e-(h3+p3)u-(hj+pj)^] 
X2+  =  exp  ZTTi  [(h2+P2)4+(h3+p3),l+(hj+Pj)4] 
X_24'  =  exp  2^i  [(h3+P3)|+(h^+Pj)ri+(h2+p3)4] 
X3+  =  exp  Eiri  [-(hj+Pj)^(h2+p3)u-(h3+p3)C] 
Yq+  =  exp  2iTi  [(hj+Pj)g+(h3+p3)Ti+(h3+p2)C] 
Y^ii*  =  exp  2iri  ['(h3+P3)4-(h2+P3)Ti-(hj+Pj)t] 
Y_j>(-=  exp  2^i  [-(h2+P2)^(hj+Pj)p-(h3+p3)^] 
Y2+  =  exp  2^i  [(h2+P2)e+(h3+pj)p+(h3+p3);] 
Y.2+=  exp  2iri  [(h3+P3)e+(h2+P2)n+(hj+Pj)C  ] 
Y3+  =  exp  2vi  [-(hj+Pj)|-(h3+p3)p-(h2+P2)& ] 
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^3d 


Table  4  (Continued) 


(RTm)  and  (RTc) 


SPACE  GROUP  cl,,  dZ,  D®,  AND  d!^.  37. 

3i  3  3a  3d 

Table  5 

Matrix  elements  and  characters  for  propagation  along  A,  z  axis.  Here 
Pi  =  P2  =  P3-  p's  equal  zero.  At  Z,  =  1/2.  Abbreviation 

u  stands  for 

d,  (Rl)  and  (R32) 


1 

X 

'  0 

2 

-2 

"^1 

1 

1 

1 

1 

u> 

2 

OJ 

^3 

1 

2 

<jj 

CO 

—  3d  D  (RTc)  The  Abbreviation  a  stands  for  e^^^^l 

For  the  group  (R3m)  a  is  to  be  replaced  by  unity. 


^0 

X3 

X-2 

^0 

^2 

Y-2 

1 

1 

1 

a 

a 

a 

^2 

1 

1 

1 _ 

1 

-a 

-a 

-a 

(Aj)ji 

I 

1 

'2 

1 

1 

"2 

a 

a 

’2 

a 

2 

^■^3^21 

n 

f  sir 

^y3 

A 

asT^ 

a»VT 

u 

-T 

■"T 

■  “Z 

"T 

<^3>12 

0 

sjT 

*T 

NfT 

2 

0 

ai\ir 

T 

o-sir 

~T 

1 

1 

~r 

1 

-a 

a 

7 

a 

2 

X  (A3) 

2 

1 

-1 

-1 

0 

0 

0 
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SPACE  GROUPS  Cjj,  D3,  Dj^.  AND 
Table  6 

Matrix  elements  and  characters  for  point  Z,  p,  =  p,  =  p,  =  1/2.  For  C?.  (RT), 

*7  5  *3^  A  «  X 

{R32),  and  (RTm),  the  symmetry  is  identical  with  r,  and  the  same 

table  can  be  used,  replacing  the  symbols  r  for  the  representations  by  Z.  The 

abbreviation  w  stands  for  e^^^^ 

d!^  .  (RTc) 


X, 

X  , 

X  0 

X., 

Y, 

Y  , 

Y^ 

Y 

1 

2 

-2 

3 

0 

1 

-1 

2 

-2 

3 

(Zj)u 

1 

0 

0 

1 

1 

0 

1 

0 

0 

1 

1 

0 

<^1)21 

0 

1 

1 

0 

0 

1 

0 

1 

1 

0 

0 

1 

<^1^12 

0 

1 

1 

0 

0 

1 

0 

-1 

-1 

0 

0 

-1 

^^1^22 

1 

0 

0 

1 

1 

0 

-1 

0 

0 

-1 

-1 

0 

X(Zi) 

2 

0 

0 

2 

2 

0 

■  0 

0 

0 

0 

0 

0 

1 

CO 

2 

CO 

2 

CO 

CO 

1 

0 

0 

0 

0 

0 

0 

^^2^21 

0 

0 

0 

0 

0 

0 

1 

-CO 

2 

-CO 

2 

CO 

CO 

-1 

<^2^2 

0 

0 

0 

0 

0 

0 

-1 

2 

-CO 

-CO 

-CO 

2 

-CO 

-1 

^^2^22 

1 

2 

-CO 

-CO 

CO 

2 

CO 

'1 

0 

0 

0 

0 

0 

0 

X(Z2> 

2 

iN/r 

-islT 

-1 

-1 

0 

0 

0 

0 

0 

0 

0 

<23)11 

1 

2 

CO 

CO 

CO 

2 

(0 

: - j 

1 

0 

0 

0 

0 

0 

0 

<23)21 

0 

0 

0 

0 

0 

0 

1 

2 

-CO 

-CO 

CO 

2 

CO 

-1 

<^3)12 

0 

0 

0 

0 

0 

0 

-1 

-CO 

2 

-CO 

2 

-CO 

-CO 

-1 

<23)22 

1 

-CO 

2 

-CO 

2 

CO 

CO 

-1 

0 

0 

0 

0 

0 

0 

XiZj) 

2 

-iN^ 

-JT 

-1 

-1 

0 

0 

0 

0 

0 

0 

0 
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SPACE  GROUPS  C^..  d!.  d!,.  AND  d!^  . 

il  3  3d  3d 

Table  7 

Matrix  elements  and  characters  for  propagation  along  A,  y  axis,  and  A', 

7 

direction  parallel  to  y  axis  passing  through  A.  For  the  group  C^.  (R3")  these 
are  not  symmetry  directions.  For  A,  p^  =0,  p^  =  -p^.  For  A',  p^  =  1/2, 
P2  =  Point  A,  Pj  =  1/2,  =  0. 


(R32),  (R7m),  and  (R76) 


Table  8 

Matrix  elements  for  propagation  along  a  direction  in  the  xz  plane.  Here  P2=p3- 
For  r,  Pj  =  P2  =  P3  =  Pj  =  ^2,  P2  =  p^  =  0.  For  B,  p^  =  0, 

p^  =  p^  =  -l/2.  For  (RT)  and  (R32)  the  general  propagation  in  the  xz 
plane  is  not  a  symmetry  direction. 

(RJm)  and  ,  (RTc) 

Even  1  1 

Odd  1  -1 
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SPACE  GROUPS  C^.,  d!,  AND  d!^^ 

3i  3  3d  3d 

Table  9 

Matrix  elements  and  characters  for  the  point  A.  Here  =  l/2,  =  0. 


"3i 


(RT) 


^0 

X3 

1 

1 

1 

-1 

dJ  (R32) 

^0 

^3 

1 

1 

^2 

1 

-1 

(RTm) 

^0 

^3 

^0 

^3 

1 

1 

1 

1 

^2 

1 

1 

-1 

-1 

^3 

1 

-1 

-1 

1 

1 

-1 

1 

-1 

Dj^  (RTc) 

^0 

^3 

^0 

Y3 

1 

0 

0 

1 

<^1)21 

0 

1 

1 

0 

<Ai)i2 

0 

1 

-1 

0 

*■^1*22 

1 

0 

0 

-1 

X(Aj)  2 


0 


0 


0 
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Table  10 

Matrix  elements  and  characters  for  the  point  B.  Here  =  0,  =  -  l/  2. 

For  (Rl*)  and  (R32)  the  symmetry  situation  is  the  same  as  at  the  point 
A.  For  (RTm)  and  (RTc)  ,  we  have  the  same  situation  as  for  (RTm) 
at  point  A.  The  matrix  elements  are  then  as  follows  for  these  two  cases: 


^0 

^3 

^0 

^3 

^1  ' 

1 

1 

1 

1 

-1 

-1 

-1 

-1 

1 

®4  ' 

-1 

1 

-1 

Table 

11 

Compatibility  Relations 

C3J  (R7) 

Pi.Zi 

^4*24 

^5-25 

^1 

^3 

^2 

^2 

^3 

D3(R32) 

r^.Zj 

^2-^2 

^3’  ^3 

^2 

^3 

^'l 

^2 

Djd  (R5m) 

r’l-Zi 

1^2’ ^2 

^3-^3 

^4-24 

r^.z^ 

^6*26 

^1 

^2 

^1 

^2 

A3 

A3 

^1 

^2 

^2 

^1 

^l'-“2 
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SPACE  GROUPS  C^,.  dI,  D®. ,  AND  d!^ 
31  3  3a  3( 

Table  11  (Continued) 


D^^(RTm)  (continued) 


A'  A'  A’ 

12  12 


D3J  (R7c) 


pl 

^2 

^3 

^4 

^5 

^6 

^2 

Aj 

^2 

^3 

A3 

^2 

^2 

^1’  ^2 

Reference 

1.  G.  F.  Koster,  Space  Groups  and  their  Representations,  Solid  State 

Physics,  Vol.  1,  p.  173,  Academic  Press  Inc.,  New  York,  1957. 
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BENZENE 


Jules  Moskowitz 


The  coafiguration  interaction  calculation  based  on  the  *'core  potential” 
approximation  discussed  in  the  past  progress  report  has  been  slightly  extended. 
The  three  center  penetration  integrals  have  been  accurately  computed  rather 
than  calculated  by  means  of  the  Mulliken  approximation.  As  anticipated,  the 
basic  results  have  not  changed  significantly.  The  more  one  improves  the  it 
electron  wave  function  the  worse  the  agreement  with  the  experimental  results. 

3 

The  state,  for  example,  is  computed  at  2.  44  ev  while  the  experimental 

value  is  3. 66  ev.  •  The  "core  potential”  approach  thus  seems  to  contain  a  funda¬ 
mental  defect. 

A  possible  remedy  would  be  the  inclusion  of  the  sigma  electrons 
explicitly.  This  is  fundamentally  a  data  processing  problem.  A  set  of  pro¬ 
grams  to  perform  the  necessary  calculations  is  in  the  process  of  being  written. 
A  self-consistent-field  calculation  on  the  benzene  ground  state  would  be  the 
initial  step  in  this  approach., 
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AROMATIC  AND  RELATED  MOLECULAR  CALCULATIONS 


B.  T.  Sutcliffe 


As  outlined  in  a  previous  QPR  (No.  47,  January  15,  1963,  p.  Il6), 
calculations  are  being  attempted  on  some  molecular  systems  using  the  method 
of  "Generalized  Product"  functions^ It  is  assumed  that  the  "group"  functions 
which  go  to  make  up  a  particular  Generalized  Product  function  are  orthogonal 
in  the  "strong"  sense  of  ref.  and  are  each  composed  of  sets  of  one  electron 
functions. 

Programs  have  been  written  which,  starting  with  integrals  over  an 
arbitrary  basis  of  one  electron  functions,  produce  integrals  over  a  suitably 
orthogonal  basis.  The  minimum  number  of  required  integrals  are  produced 
and  they  are  arranged  in  such  a  way  that  the  elements  of  the  effecfive  Hamil¬ 
tonian  (see^^^  eq.  25,  p.  363)  for  each  group  in  turn  may  be  made  up  quickly 
and  easily. 

It  is  hoped  shortly,  to  perform  a  pilot  calculation  on  ethylene  using 
Gaussian  functions  to  obtain  the  basic  integrals. 

Reference 

1.  R.  McWeeny,  Rev.  Mod.  Phys.  32,  335,  (I960). 
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THE  POLYATOM  SYSTEM 


M.  C.  Harrison 


The  POLYATOM  system  of  programs  in  its  present  preliminary  form 
is  capable  of  evaluating  LCAO-SCF  wavefunctions  for  a  closed  shell  molecular 
system,  expressed  in  terms  of  not  more  than  50  s  and  p-type  Gaussian 
functions  centered  on  not  more  than  20  nuclei.  The  various  programs  com¬ 
prising  the  system  communicate  with  each  other  via  magnetic  tape,  on  which 
data  is  stored  in  a  specified  format.  It  is  convenient  to  crystalize  this  format 
at  an  early  stage  of  development,  so  that  programs  using  the  same  data  may  be 
written  in  parallel.  Accordingly,  the  content  and  format  of  two  of  the  mag¬ 
netic  tapes  if  given  in  detail  below.  It  should  be  noted  that  the  file-handling 
routines  described  briefly  in  CCL  Technical  Note  19  are  used  in  referring  to 
magnetic  tape.  This  means  essentially  that  each  file  of  information  is  given 
a  BCD  name,  and  can  be  located  on  the  tape  when  this  name  is  known,  what¬ 
ever  order  the  files  may  be  in. 

(1)  Tape  KMLIST 

This  tape  is  a  logical  unit  12,  and  contains  3  files,  named  K-LIST, 
V-LIST,  and  M-LIST.  Given  an  ordered  set  of  one -electron  functions 
^2'  *  *  ■  certain  symmetry  properties,  these  files  contain  labels  of 

non-zero  integrals  of  the  types  (t).  (rj.  |v|r|^),  and  where 

K  is  any  totally  summetric  one-electron  operator  such  as  unity  or  -l/2  V 
is  any  operator  having  the  symmetry  of  the  nuclei  such  as  S  —  ,  and  M  is  any 

1  ^ 

totally  symmetric  two-electron  operator  such  as  —  .  The  indices  in  each  label 
are  given  in  standard  order;  that  is,  with  i^j  in  tSe  one-electron  cases  and 
with  i  >  j,  k  >  i,  and  (ij)  >  (ki)  in  the  two-electron  case,  where  (ij)  =  i/2  (i  +  1)  +  j. 
The  labels  of  integrals  which  are  identically  equal  within  a  minus  sign,  because 
of  the  symmetry  properties  of  the  tiSjare  given  consecutively. 
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Each  of  the  three  files  is  divided  into  Fortran  logical  binary  records, 
with  the  first  word  in  each  record  being  a  Fortran  integer  not  greater  than  500 
giving  the  number  of  labels  in  the  record,  the  second  a  Fortran  integer  which 
is  1  if  the  record  is  the  last  in  the  file,  and  zero  otherwise,  and  the  remainder 
the  labels,  one  in  each  word.  Each  word  is  divided  into  6  fields,  comprising 
bits  S“5,  6-11,  12-17,  18-23,  24-29,  30-35,  For  the  one  electron  integrals, 
fields  1,  2,  and  6  are  used  to  store  i,  j,  and  a  tag  respectively.  For  the  two- 
electron  integrals,  fields  3  and  4  are  used  to  store  k  and  i.  The  value  of  the 
tag  is  zero  if  the  integral  is  not  identically  equal  to  any  preceding  integral,  one 
if  it  is  equal  to  the  last  integral  with  a  zero  tag,  and  two  if  it  is  equal  to  minus 
the  last  integral  with  a  zero  tag. 

The  ordering  of  the  labels  is  such  that  the  integrals  which  are  equal 
within  a  minus  sign  are  in  standard  order  (that  is,  in  increasing  order  of 
((ij)  (ki) )  ),  and  the  integrals  with  zero  tag  are  also  in  standard  order. 

(2)  Tape  KMINTS 

This  tape  is  a  logical  unit  11,  and  contains  4  files,  named  G-INTS, 
T-INTS,  V-INTS,  and  M-INTS,  These  contain  respectively  the  values  of  the 
overlap,  kinetic  energy,  potential  energy,  and  electronic  repulsion  integrals. 
Each  file  is  divided  into  Fortran  logical  binary  records,  with  the  first  two 
words  in  each  file  having  the  same  meaning  as  in  (1)  above.  The  remainder 
of  the  record  consists  of  the  labels  as  before,  but  with  each  label  followed  by 
the  value  of  the  integral  as  a  floating-point  word.  No  assumptions  should  be 
made  about  field  6  of  the  label,  however.  It  should  not  be  assumed  that  there 
is  only  one  integral  for  given  values  of  the  indices.  Thus  file  V-INTS,  for 
instance,  may  contain  several  entries  for  integral  (i  |v  jj) ,  corresponding  to 
different  centres,  which  should  be  added  together. 

Note  that  the  present  format  can  only  be  used  for  up  to  63  one -electron 
functions.  Should  this  number  be  insufficient  the  tape  format  can  easily  be 
changed.  However,  preliminary  tests  on  the  709  suggest  that  the  computation 
time  necessary  for  a  larger  basis  set  will  be  excessive. 


MOLECULAR  INTEGRALS  OVER  SLATER  ORBITALS 


D.  E.  Ellis 


I.  One  and  Two  Center  Integrals 

A  set  of  computer  routines  has  been  prepared  by  Mr.  C.  A.  Christy 
and  the  author  for  the  computation  of  all  one  and  two  center  integrals  arising 
in  energy  calculations  using  a  basis  of  Slater  orbitals.  The  one  center  inte¬ 
grals  are  computed  from  analytic  formulae,  and  two  center  integrals  are 

(2) 

evaluated  by  the  Barnett- Coulson  method'  All  angular  integrations  are 

carried  out  analytically.  The  radial  integrations  are  performed  by  making 

(3) 

use  of  the  analytic  series  expansions  for  molecular  zeta  functions'  ’ 
except  for  the  exchange  integrals,  <  AB  |  AB  >.  The  <  AB  |  AB  >  integrals 
appear  as  an  infinite  series  of  radial  integrals  multiplied  by  angular 
coefficients,  and  these  radial  integrals  are  computed  by  means  of  Gauss- 
Legendre  and  Gauss-Laguerre  quadratures^^^. 

Testing  of  the  routines  has  been  completed  for  combinations  of  Is, 

2s,  2p,  and  in  some  cases,  3d,  orbitals.  Results  have  been  compared  with 

(4) 

the  DIATOM  program  written  by  Switendick  and  Corbato'  and  with  several 
other  integrals  programs.  Further  testing  is  planned. 

Some  comparison  may  be  made  with  the  DIATOM  program,  which 
has  been  in  frequent  use.  The  angular  integrations  proceed  in  much  the  same 
manner  in  both  programs.  DIATOM  makes  use  of  equal  interval  rule  numeri¬ 
cal  integrations  for  all  radial  functions,  including  one  center  integrals.  The 
analytic  and  Gauss  integrations  used  in  the  routines  being  described  permit 
higher  accuracy  and  greater  speed.  The  DIATOM  program  computes  all 
non-zero  integrals  for  a  given  basis  set,  and  allows  the  basis  set  to  be  chosen 
as  a  linear  combination  of  Slater  orbitals.  Sometimes  only  a  few  integrals, 
or  integrals  of  a  particular  type  are  desired  and  it  would  be  convenient  to  be 
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able  to  compute  only  those  integrals.  A  very  simple  main  program  making 
use  of  the  set  of  routines  is  being  written  for  this  purpose. 

In  the  computation  of  a  large  number  of  integrals  required,  for 
example,  in  a  self-consistent  molecular  calculation  it  becomes  possible  to 
decrease  the  computation  time  per  integral  very  considerably  by  saving 
auxilliary  functions  from  one  integral  to  the  next.  As  an  example,  exchange 
integrals  requiring  1.  5  minutes  computation  singly  have  been  produced  at  an 
average  rate  of  about  10  sec.  apiece  when  a  number  are  computed  at  once. 
Some  function  saving  has  been  incorporated  into  the  set  of  routines,  and  a 
main  program  is  planned  which  will  set  up  an  efficient  calculation  of  integrals 
over  an  entire  basis. 

IL  Three  Center  Integrals 

A  computer  routine  has  been  written  and  partially  tested  for  the 
evaluation  of  the  three  center  nuclear  repulsion  integrals,  <  AbI^^  >.  This 
routine,  together  with  the  routines  described  above,  allow  the  computation 
of  all  integrals  appearing  in  the  one-eTectron  Hamiltonian.  Formulae  have 
been  derived  for  the  two  electron  coulomb,  <  AA  |  BC  >  and  exchange 
<  AB  I  AC>  integrals,  and  are  being  programmed.  An  integrals  package  is 
planned  which  will  combine  all  routines  under  control  of  a  main  program  for 
most  efficient  calculation. 

III.  Four  Center  Integrals 

Elsewhere  in  this  report  Mr.  J.  P.  Wright  describes  a  program 
which  he  is  writing  to  evaluate  four  center  electron  repulsion  integrals  be¬ 
tween  combinations  of  s,  p,  and  d  orbitals.  We  plan  to  co-ordinate  the 
programs  described  above  with  this. 

IV.  Formulae 

We  want  to  sketch  the  derivation  and  give  the  formulae  for  a  few  of 
the  integrals,  but  for  the  sake  of  brevity  some  steps  and  definitions  will  be 
omitted. 

Define  the  normalized  Slater  orbital  as 

.  n-1  -Kr.  .m_,m/  v  o  /x\ 

<j>  =  N  r  e  (-) 

by  Hobson's  definition,  m  >  0 
S^((t))  =  cos  (m<t>) 
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arid  use  the  following  expansions  and  definitions 


m-1  -/cr. 


i  -dHi  ,  v  /X  \  V  /  vi+m.i+mv  i  xrycr  ,  ,  . 

<%)  =Z  <i+m)  Pi 


(x-axes  antiparallel) 


00 

J_  =  V  6.(r  ,  r  )  P^{ 

''12  fZ  j  ^1  ^2  ^ 


These  relations  prove  sufficient  for  the  derivation  of  all  one  and  two 
center  formulae,  and  the  three  center  nuclear  repulsion  integral.  It  is  assumed 
that  the  coordinate  system  on  center  A  is  right  handed,  and  that  system  B  is 
left  handed,  with  z-axes  antiparallel.  To  allow  parallel  co-ordinate  axes  in  an 
absolute  frame,  and  to  compute  the  remaining  three  and  four  center  integrals 
requires  the  introduction  of  the  rotation  matrices  for  solid  harmonics.  This 
we  delay  until  a  future  report. 

The  method  of  derivation  is  to  expand  all  functions  onto  center  A,  to 
couple  the  resulting  Legendre  functions,  to  expand  products  of  trigonometric 
functions,  and  to  perform  the  remaining  angular  integrals  analytically.  One 
important  result  is  an  expression  for  the  potential  integral 


<Ab|Z->=  ^.(A,  1)4>2(B.  1)7^  dv 

^12-'^  ^  ^12  ^ 


-  I  I  JFiI 

J=|mj-m2l  l2= 


IJTlZ  <->  <i^m 


l2=m2 


^  (2X2+1  ){ang||z 
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r  I  r  12  2'  1-  2  -  '1  2' 

where  |aag}-=|(-)  S 


m^+m^  I  2 


^2*  ^I'  ^2  ’^2  *  “^2 


+  s 


Im^-m^  I 


(the  F's  are  a  sum.  of  products  of  two  O^’s) 


and 


iY}=^e  '"2  t  ^  ^ 


To  obtain  the  three  center  nuclear  attraction  integral,  we  merely  set  r^  =  R^, 

4>=#=0fx  =M'  -  The  two  center  exchange  integral  is  obtained  ai 

a  c  ^7  ^ 


follows: 

<  AB  AB  >  =  \  <  Afi.  —  >  AB  dv- 
J  rj2  ^-2 


and  the  result  is 


< Ab|aB  >  =  6  I 

'”'r‘^2 


Y  V"'4/V’"4\  ^4'^4V  i2\2+^){Z^,  +  l)  f_  1  r  1 


4  4 


4-^4 

finite 


where 


4=1 


KZ  ^3+14+1-  «4 

zx  f-=  \  e  s  ^  .  X  (1  •—  s, 

0  n4-^4->^4  ^2 


"4«b)  I 


Kj  n.+i-  +  l 
e  ^  t  ^  ^ 


and  ^  ang  is  a  sum  of  four  terms,  each  of  which  is  a  product  of  two  F's  and 
a  factorial  ratio. 
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J.  P.  Wright 


Testing  has  begun  of  a  program  written  to  calculate  four-center 
molecular  integrals  for  any  combination  of  s,  p,  and  d  functions  with  un¬ 
restricted  principal  quantum  numbers.  The  zeta  function  method  of  Barnett 
and  Coulson  has  been  used.  The  formula  for  the  general  case  has  been  derived 
using  a  notation  similar  to  that  given  by  D.  E.  Ellis  in  a  report  elsewhere  in 
this  QPR,  where  he  discusses  his  work  on  three -center  molecular  integrals. 
The  formula  will  be  given  in  a  future  report,  after  testing  has  been  com¬ 
pleted. 
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EXPANSION  THEOREMS  FOR  SOLID  SPHERICAL  HARMONICS 


J.  P.  Dahl 
M.  P.  Barnett 


Some  expansion  theorems  for  solid  spherical  harmonics  have  been 
derived.  It  is  hopeful  that  these  theorems  will  prove  to  be  useful  in  the  calcu¬ 
lation  of  three  and  four-center  molecular  integrals. 

Consider  the  two  coordinate  systems  a  and  b  (Figure  1)  with  parallel 
axes.  The  location  of  b  relative  to  a  is  specified  by  giving  the  vector 
R  =  (X,  Y,  Z)  =  (R,0,^,  where  the  rectangular  and  spherical  coordinates  of  R 
expressed  in  a  have  been  denoted  (X,  Y,  Z)  and  (R,  9,^)  respectively.  With  a 
similar  notation  a  field  point  P  may  be  specified  by  giving  its  coordinates  in 
b.  i.e.  =  r  =  The 

relation 

(x,  ,y,  ,z  )  =  (x  -X,  y  -Y,  z  -Z)  (1) 

b  •'b  b  a  ‘^a  a  ' 


obviously  holds. 


p 


Figure  1 
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The  following  expansions  are  then  valid: 


rN  pM  (eos 


e-iM4.b 

cos<^<t>b) 


N  n 

y  y  (-)^'"  r'^  P“(cos  e  )  p|^-’”(cos  0) 

LiLi  ^  ^^+^1  /a  n  '  a'  N-n  '  ' 

n=0  m= 

n-N+M 


im9  i(M“m)<P 

^  e  (2) 

-im<|)  -i{M-m)# 

^  e  (3) 


cos[m4)g^  +  (M-m)#]  (4) 

sin[m4>g^  +  (M-m)#]  (5) 

The  associated  Legendre  functions  are  those  defined  by  Hobson^ 

,m 

P^(cos  ^  =  (-)^  sin”^0 - —  P  {cos  m  >  0  (6) 

^  d{cos  ^ 


„-m ,  ^ 

P^  (cos  6)  =  (-) 


(n-m)! 

(n+rn)T 


P“'  (cos  ^ 
n 


(7) 


The  relation  (7)  holds  for  all  m  values.  Further  it  follows  from  (6) 
that  p|^(cos  6)  =  0  for  |m  |  >  n. 

The  coefficients  occurring  in  (2)  *  (5)  are  the  binomial 

coefficients,  i.  e.  : 

/P\  -  P! 

“  g-np-g)l 

It  suffices  to  demonstrate  the  validity  of  (2),  since  (3)  is  merely  the 
complex  conjugate  of  (2).  Further  (4)  and  (5)  are  obtained  from  equations  (2) 
and  (3)  simply  by  adding  and  subtracting  these  two  equations. 

From  the  defining  equation  (6)  it  follows  that 

rJi  ,  ^  .n  (2n)!  .  n^ 

P  (cos  6)  =  (1)  -  sin  G 

n  I 

2  n  ! 
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^  inA  /  »n  (2n)|  .  .  • 
r  P  (cos  6)  e  ^  =  (-)  ■  '  (x+iy) 


-n  , 
2  n  ! 


(8) 


Coasidering  Eq.  (1)  we  then  obtain 


N 


,.,N  ,^y  ,.,N-n,N„  ,  ,.,x«Y,«-n 

n  a  a 

n=:0 


n=0 


N 


in(})  i(N-n)# 

e  ^  e 


(9) 


Equation  (9)  is  obviously  identical  with  equation  (2)  in  the  case  M-N. 

Equation  (2)  has  thus  been  shown  to  hold  for  the  maximum  value  which  M  can 

take.  From  the  process  of  induction  it  follows  that  the  equation  is  valid  in 

general,  if  validity  for  M  implies  validity  for  the  value  M-1 . 

Suppose  therefore  that  (2)  holds  for  the  value  M.  We  then  consider 

simultaneous,  infinitesimal  rotations  of  the  coordinate  systems  a  and  b 

through  the  same  angle  about  parallel  axes.  The  changes  in  the  functions  in- 

volved  are  then  found  by  applying  the  operators  of  angular  momentum,  since 

these  operators  are  proportional  to  the  operators  for  small  rotations.  The 

fact  that  r^,  r^^  and  R  are  all  rotated  through  the  same  angle  then  implies  that 

operating  with  an  angular  momentum  operator  on  the  left  hand  side  of 

Eq.  (2)  is  equivalent  to  operating  with  the  operator  i  +  i  on  the  right  hand  side. 

a  ^ 

Here  I  and  are  the  same  functions  of  (0,  #)  and  3is 

From  the  well  known  properties  of  the  operator  -  ii^  it  follows 

that 

(f  -if„)  P^(cos  G)  =  (n+m)(n“m+l)P^  ^(cos  ^  (10) 

X  y  n  n 

Operating  on  the  left  hand  side  of  Eq.  (2)  with  the  operator  i  -  ii 

and  on  the  right  hand  side  with  the  operator  (i  -ii^)  +  (i  -  ii  )  we  obtain: 
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(N+M)(N-M+l)rj^P^‘Ncos 


I  I  (-)”•“  O'™''"-*" 

a=0  •m= 

a-N+M 


p;^-l(cos  0^)p;y.7(cos  0)e 


i(M-m)# 

e 


Y  (^+m) 

n=0  m= 

m-N+M 


P 


m 

n 


(cos  0^)  P 


M-m-l 

N-n 


im(t)a 

(cos  0)  e 


i(M-m-l)# 

e 


Writing  m  for  m-1  in  the  first  summation  on  the  right  hand  side  it  is 
readily  shown  that  this  equation  may  be  rewritten  as 

NpM-l,  .  i(M-l)(t>^ 


(-) 


,N-n  /N+M-l\  ^n  „N-n 


\  n+m 


r“  R" 
a 


n=0  m= 

n-N+M-1 


Tv/r  1  i(M-m-l)# 

p"^(cos  0  )  pM-m-1  g 

n  '  a  N-n  ' 


which  is  just  Eq.  (2)  for  the  value  M-1 . 

The  general  validity  of  the  equations  (2)  -  (5)  has  thus  been  demon¬ 
strated. 

After  completion  of  the  present  work  it  has  come  to  our  attention  that 
an  equation  similar  to  Eq.  (2)  has  been  derived  by  M.  E.  Rose  (J.  Math,  and 
Phys.  37,  215  (1958)),  who  applied  the  algebra  of  irreducible  tensors.  The 
equation  appeared  as  a  subordinate  result  in  a  study  of  the  electrostatic  inter¬ 
action  of  two  charge  distributions.  Rose  presented  his  equation  in  the  form: 


N 


I  I  1-1'''^ ,R)  Ir- I 


n=0  m=-ii  [(2n+l)  !  (2N-2n+l)!  ] 


1/2 


where  is  a  normalized  solid  spherical  harmonic; 
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r«P;^(cos  a 

and  the  C-coefficients  are  vector  addition  coefficients.  By  substituting  the 
actual  form  of  these  coefficients  into  Rose’s  equation  the  equivalence  with 
equation  (2)  could  be  demonstrated. 
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MONTE  CARLO  METHOD  FOR  MOLECULAR  INTEGRALS 


C.  W,  Nielson 
B.  J.  Woznick 


The  Monte  Carlo  method  has  been  widely  used  for  solving  problems  in 
neutron  kinetics  and  has  sporadically  been  used  for  the  evaluation  of  multi¬ 
dimensional  integrals  for  other  purposes 2,  3,  4)^  investigating  the 

possible  usefulness  of  this  method  for  the  evaluation  of  multicenter  molecular 
integrals. 

The  advantages  of  this  method  are  its  extreme  simplicity,  the  fact 
that  it  is  not  limited  to  any  special  type  of  orbital  functions,  and  the  fact  that 
it  provides  automatically  an  estimate  of  the  error.  However,  Monte  Carlo  cal¬ 
culations  are  almost  certain  to  take  much  longer  for  comparable  accuracy  than 
methods  which  make  special  use  of  the  analytic  properties  of  the  functions  in¬ 
volved,  such  as  the  zeta-function  expansion  method  or  the  Gaussian  transform 
method.  Nevertheless,  it  seems  that  the  Monte  Carlo  approach  has  real  value 
for  the  evaluation  of  integrals  involving  unusual  functions,  perhaps  numerical, 
and  also  for  the  checking  of  other  methods  as  they  are  extended  into  new  do¬ 
mains.  In  addition,  the  simplicity  of  the  method  makes  it  feasible  to  implement 
it  on  very  small  computers. 

The  basis  of  the  method  may  be  seen  by  considering  the  definition  of 
an  Integral  as 

f  f(x)dx  =  ^  f(xj  AXk  (1) 

^  Axk-0  k=l 

where  Ax^  =  is  some  point  in  this  interval.  Then,  if  the 

intevals  are  chosen  to  be  equal  and  the  limits  taken  as  0  and  1, 


(2) 
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■  ^ 

wh'er^e  the  are  distributed  uniformly  in  the  interval  0-1.  Ordinary  techniques 
of  integration  achieve  this  uniform  distribution  by  chopping  the  interval  into  n 
equal  subintervals.  The  essence  of  the  Monte  Carlo  method  is  to  choose  these 
points  at  random;  in  this  context,  "at  ramdom"  simply  means  that  we  have  a 
set  of  points  which  are  uniformly  distributed  and  statistically  independent. 

There  is,  of  course,  no  value  in  doing  this  for  one -dimensional  inte¬ 
grals.  However,  theorems  of  probability  theory  have  been  proven  which  show 
that  the  probable  error  in  the  integral  goes  down  as  l/^yN  regardless  of  the 
dimensionality.  In  contrast,  straightforward  numerical  integration  requires 
the  number  of  points  to  increase  geometrically  with  the  dimension.  For  a 
sixfold  integration  such  as  occurs  in  molecular  integrals,  straightforward 
numerical  integration  seems  out  of  the  question,  whereas  Monte  Carlo  integra¬ 
tion  does  not. 

Another  important  tool  in  making  the  technique  practicable  is  the  use 
of  a  distribution  function,  or  weighting  factor.  Assume  that  the  integrand  can 
be  written  as  the  product  of  two  factors 


f(x)  g(x)  dx 


(3) 


where 


I 


g(x)  dx  =  1 


(4) 


and  g(x)  is  everywhere  non-negative.  Then  we  can  take  this  function  to  be  a 

❖ 

probability  density  and  select  the  points  according  to  this  distribution. 
This  having  been  done, 


f(x)  g(x)  dx  = 


Urn 

N-*-oo 


Li 


(5) 


Qualitatively,  this  procedure  can  be  looked  upon  as  dividing  out  an  exactly 
integrable  part  of  the  integrand  leaving  a  more  smoothly  varying  function  to 
integrate  numerically.  If  the  distribution  function  is  chosen  appropriately,  the 
numerical  integration  points  are  taken  only  where  the  integrand  is  significant, 
hence  the  technique  is  sometimes  called  importance  sampling. 

We  have  used  these  methods  to  calculate  four-center  molecular  inte¬ 
grals  using  IS  orbitals  on  all  centers.  The  distribution  function  used  was  a 
product  of  two  three-dimensional  factors,  one  for  each  of  the  two  electronic 
coordinates  systems  of  the  molecular  integral.  Each  factor  was  a  spherically 
symmetric  function  centered  somewhere  in  the  molecule.  The  radial  part  of 
the  function  was  a  constant  inside  a  radius  characteristic  of  the  molecule  and 


of  the  form  e 


outside  this  radius  where  t  was  chosen  to  match  the 


r 
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assymptotic  form  of  the  orbitals.  This  function  is  suitable  for  a  distribution 
function  because  it  can  be  integrated  exactly  and  thus  normalized.  Random 
numbers  uniform  in  the  unit  hypercube  are  generated  on  the  computer  by 
standard  algorithms;  these  are  more  properly  called  pseudo-ramdom  numbers, 
since  they  are  in  fact  generated  by  an  exact  algorithm.  They  do,  however,  share 
most  of  the  mathematical  properties  of  random  numbers.  The  uniform  distri¬ 
bution  was  transformed  by  a  "stretching”  process  subject  to  the  requirement 

dx  =  r  f(r')dT  (6) 

0  *^0 

where  p  is  the  radius  in  the  unit  sphere  and  r  is  the  radius  in  the  final  distri¬ 
bution.  (Points  falling  outside  the  unit  spheres  are  discarded.) 

As  an  example  of  the  results,  a  typical  integral  was  obtained  using 
5000  random  points  as  .35  ±  .01  whereas  the  zeta-function  programs  availa¬ 
ble  give  in  a  comparable  time  0. 36735,  a  value  believed  to  be  correct  to  all  of 
the  figures  given.  The  Monte -Carlo  computer  program  required  for  this 
calculation  is  extraordinarily  short  compared  to  the  zeta-function  program 
system.  Also,  negligible  amounts  of  temporary  storage  are  needed.  Thus  it 
is  demonstrated  that  integrals  of  limited  accuracy  can  be  obtained  very  simply 
by  this  method.  However,  in  view  of  the  fact  that  the  Monte -Carlo  error  would 
decline  only  as  1/ n/N,  accuracy  comparable  to  the  zeta-function  method  is 
practically  impossible. 

(5) 

Recently  Haselgrove'  '  has  suggested  using  so  called  diophantine 
numbers  in  a  fashion  similar  to  the  random  numbers  in  the  Monte -Carlo 
method.  These  numbers  also  tend  towards  filling  space  uniformly,  but  with 
less  fluctuation  from  uniformity  than  random  numbers.  His  claim  is  that  the 
error  in  an  integral  in  this  case  goes  down  as  l/ N  rather  than  l/ n/N  as  in  the 
Monte  Carlo  case.  We  are  now  investigating  the  application  of  these  dio¬ 
phantine  numbers  to  the  evaluation  of  molecular  integrals. 
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THE  QUADRATIC  JAHN-TELLER  EFFECT 


Don  E.  Ellis 
Robert  Englman 


The  ground -state  energy  expectation  values  have  been  computed,  as 
functions  of  the  strength  of  the  linear  interaction  (k)  and  the  quadratic  coupling 

(i),  from  the  Hamiltonian  and  the  approximate  solution  given  in  an  earlier 
report^ The  questions  that  the  computation  was  designed  to  answer  were: 

(1)  What  is  the  effect  of  the  quadratic  coupling  on  the  stabilization  of  the 
system  (i.  e. ,  on  the  lowering  of  the  ground  state  energy  compared  to  the 
energy  of  an  uncoupled  system)? 

(2)  Is  any  analytic  expression  suggested  by  the  results? 

(3)  What  is  the  nature  of  the  splitting  of  the  lowest  pair,  which  is  degenerate 
in  the  linear  effect  (i=0)?  In  particular, 

(4)  What  is  the  balance  of  the  effects  on  the  upper  of  the  pair  of  the  upward 
tendency  of  the  splitting  and  of  the  downward  tendency  of  stabilization,  and 

(5)  Where  does  this  level  meet  the  second  harmonic  ? 

Figures  1-4  provide  the  answer.  It  is  important  to  note  that  the 
figures  are  plotted  as  though  curves  with  different  i's  coincided  in  k=0.  This 
is  not  the  case,  since  the  z6ro-point  motion  energies  are  ^  {^JTTT  +  sfl  i  )  for 
k=0.  This  slowly  varying  quantity  should  then  be  added  to  the  value  of  the 
curve  to  obtain  the  energy. 

Regarding  question  (2),  the  plotted  stabilized  energy  can  be 
represented  by  -i  k^/  (1-/)  +  0(1)  for  large  k  and  any  i.  No  simple  expression 
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was  found  for  the  upper  of  the  pair  in  this  limit.  However,  an  analytic,  exact 
expression  was  derived  by  second -order  perturbation  theory  and  valid  up  to  the 
order  of  (for  small  k)  and  any  i  for  both  of  the  pair.  The  expression  is  in¬ 
volved.  For  i  small  it  reduces  to  (1-i)  for  the  lower  of  the  pair  and 
-k^/  (1+i)  for  the  upper  (i>0)  . 

Figure  4  is  purported  to  answer  (5).  Our  method  applies  only  to  the 
ground -state  pair  and  cannot  give  the  level  of  the  second  harmonic.  However, 
we  have  reason  to  expect  that  the  separation  of  this  level  from  the  lowest  one 
varies  only  weakly  with  k,  i  and  can  be  taken  as  1. 

The  table  gives  numerical  results  for  comparison  with  future 
computations.  Our  computations  were  carried  out  on  the  IBM  709  by  means  of 
numerical  quadrature  of  complicated  formulae.  Characteristically,  a  thousand 
values  were  computed  on  a  rectangular  mesh,  a  pair  of  energies  was  obtained 
in  2  minutes. 

The  importance  of  this  work  may  be  seen  by  the  light  of  the  fact  that 
after  the  initial  impetus  of  the  Jahn-Teller  effect  in  1957-8,  progress  has 
slowed  down  and,  in  terms  of  exact  results,  virtually  halted,  ostensibly 
because  of  the  difficu:+i?s  in  obtaining  exact  solution.  Regarding  the  reliability 
of  our  results,  based  as  they  are  on  approximate  formulae,  we  are  confident 
that  they  reproduce  faithfully  any  qualitative  aspect  and  suffice,  at  this  stage, 
for  comparison  with  th*-"^  experiments  that  one  would  think  of. 

We  are  thankful  to  Professor  W.  Thorson  for  directing  our  attention 
to  this  work. 


Table  1 

Numerical  Level  Heights  for  i  =  0.  5 


k  = 

0.  3 

0.  5 

1.0 

2.  0 

Upper 

-0. 009714 

-0. 227310 

+0.  71908 

+16. 4301 

Lower 

-0. 108763 

-0.  301881 

-1. 12986 

-4. 1443 
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<=0.9  <=0.75  t=Q5  t=0.3  t=0.l 


Figure  1 


Energy  vs.  k  (the  linear  coupling  strength)  with  i  (the  quadratic  coupling 
strength)  as  parameters.  The  drooping  curves  are  the  ground -state 
energies;  the  rising  ones  belong  to  the  upper  of  the  ground-state  pair. 


(Note  the  sliding  zero  of  the  ordinate. 
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Figure  2 

Enlargement  of  Figure  1  for  small  k 
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k  =  0.75 
k  =  0.50 


k  =  0.50 
k=0.75 

k=  1.0 


Figure  3 

Energy  vs.  i  with  k  as  parameters 
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Figure  4 


The  value  of  k  such  that  the  splitting  of  the  ground -state  pair  levels  is  1, 
plotted  vs.  i.  At  the  plotted  points  (i,  k)  the  upper  level  coincides  approxi¬ 
mately  with  the  second  harmonic  of  the  ground  state. 


